Pis'ma v ZhETF 



(N 

o 

O 

Q 

o 



6: 

cd. 

o: 

>> 
^: 

Or 



Direct numerical experiment on measuring of dispersion relation for 
gravity waves in the presence of condensate 

A. O.Korotkevich+* 1 '> 

'Department of Mathematics and Statistics, University of New Mexico, MSC01 1115, 1 University of New Mexico, Albuquerque, NM 

87131-0001, USA 

*L.D. Landau Institute for Theoretical Physics RAS, 2 Kosygin Str., Moscow, 119334, Russian Federation 

Submitted 2 December 2012 

During previous numerical experiments on isotropic turbulence of surface gravity waves we observed for- 
mation of the long wave background (condensate). It was shown (Korotkevich, Phys. Rev. Lett. 101(7), 074504 
(2008)), that presence of the condensate changes a spectrum of direct cascade, corresponding to the flux of 
energy to the small scales from pumping region (large scales). Recent experiments show that the inverse cas- 
cade spectrum is also affected by the condensate. In this case mechanism proposed as a cause for the change 
of direct cascade spectrum cannot work. But inverse cascade is directly influenced by the linear dispersion 
relation for waves, as a result direct measurement of the dispersion relation in the presence of condensate is 
necessary. We performed the measurement of this dispersion relation from the direct numerical experiment. 
The results demonstrate that in the region of inverse cascade influence of the condensate cannot be neglected. 
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Theory of wave or weak turbulence 1 applied to 
gravity waves on the surface of the fluid is the base 
for all current wave forecasting models, which are cru- 
cial for ocean cargo communications and oil and gas sea 
platforms operations. This is why it's verification is im- 
portant and urgent problem. Numerous attempts to get 
a spectrum of the direct cascade of energy from ocean 
and sea observations give results which confirm the wave 
turbulence theory [21 [3] ■ At the same time all these 
experiments were working with wind generated waves, 
which means broad spectrum pumping. As a result it 
is hard to understand where we have so called "inertial 
interval" (region of scale where we have only nonlinear 
interaction of waves, while pumping and damping influ- 
ences are negligible). Narrow frequency pumping can 
be realized in experimental wave tanks and flumes. But 
this state of the art experiments were producing very 
strange and contradictory results [H[S] , like dependence 
of the spectral slope on the amplitude of pumping. The 
direct numerical simulation of the primordial dynami- 
cal equations looks like a natural remedy for this prob- 
lem. It provides us with all possible information about 
the system, but for the cost of enormous computational 
complexity. 

One of the first attempts was pioneering work [6], 
which soon was followed by [7H9] During the last decade, 
the author together with colleagues were able to find 
answers at least to some of the open questions using 
direct numerical simulation of gravity waves. It was 
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shown, that on a discrete grid of wavenumbers (com- 
mon situation for both pseudo spectral numerical codes 
in a periodic domain and experimental wave tanks which 
are usually rectangular basins of finite size) the meso- 
scopic turbulence can take place [10l[TT]. In the recent 
works [THE!] t ne author demonstrated that formation 
of the inverse cascade, corresponding to the flux of wave 
action (analog of number of waves), inevitably leads 
to the formation of the strong long wave background, 
which we call condensate (due to similarity with Bose- 
Einstein condensation in condensed matter physics). It 
was shown that presence of the condensate changes the 
slope of the direct cascade spectrum. In these recent ex- 
periments, in spite of the short inertial interval for the 
inverse cascade, the slope significantly different from the 
predicted by the wave turbulence theory was observed. 
Recent reports |14] with long enough inertial interval 
show that the slope is indeed differs from the theoreti- 
cally predicted one. 

In the present Letter we report results of direct mea- 
surement of the dispersion relation of the surface grav- 
ity waves in the presence of condensate. As it will be 
demonstrated later, the slope of the spectrum of the 
inverse cascade directly depends on the power of the 
dispersion relation. Although we were not able to deter- 
mine the change of this power, we demonstrate that in 
the region of the inverse cascade the dispersion relation 
is strongly affected by the nonlinear interaction with the 
condensate. This distortion cannot be neglected and 
might be considered as a possible cause of the change of 
the slope of the inverse cascade spectrum. 
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As in many previous works 0[§l[T2J[T2J[Ii] we shall 
consider isotropic turbulence (no direct dependence on 
angle) of gravity waves on the 2D surface of the 3D fluid. 
The isotropic turbulence is a wonderful sandbox for the 
turbulence simulation. In this case we can use angle 
averaging of the resulting spectra in order to decrease 
it natural peakedness of harmonics, instead of averag- 
ing over different realizations which is unrealistic in this 
case, due to a long time of even one simulation. 

Here and further we shall follow notations from [T2] . 
We consider a potential flow of ideal incompressible 
fluid. System is described in terms of weakly nonlinear 
equations [I] for surface elevation rj(r, t) and velocity 
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Here dot means time-derivative, A — Laplace opera- 
tor, k is a linear integral operator {k = \J — Aj , F^ 1 is 
an inverse Fourier transform, 7^ is a dissipation rate 
(according to recent work [TB] it has to be included 
in both equations), which corresponds to viscosity on 
small scales and, if needed, " artificial" damping on large 
scales. Pf is the driving term which simulates pump- 
ing on large scales (for example, due to wind). In the 
fc-space supports of 7^ and P^ are separated by the in- 
crtial interval, where the Kolmogorov-type solution can 
be recognized. These equations were derived as a results 
of Hamiltonian expansion in terms of krj. From phys- 
ical point of view fc-operator is close to derivative, so 
we expand in powers of slope of the surface. In most of 
experimental observations average slope of the open sea 
surface \i is of the order of 0.1, so such expansion is very 
reasonable. Additional details can be found in |7ll74 - fT5] . 

In the case of statistical description of the wave field, 
Hasselmann kinetic equation [2UJ for the distribution of 
the wave action n(k,t) ~ (|a.g(i)| 2 ) is used. Here 
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are complex normal variables. For gravity waves Uk — 
\fgk. In this variables, if we have a linear wave with 
wavenumber k, it will correspond to the only excited 
harmonics ar. In other words, representation in terms 
of these normal variables means representation in terms 
of elementary excitations in the system (linear waves). 



In reality we should use pair correlator for variables af- 
ter canonical transformation which eliminates nonreso- 
nant terms in the Hamiltonian [TJ[T7], but in the case 
of gravity waves of average steepness \/(|V??| 2 ) ~ 0.1 
their relative difference is of the order of few percents. 
Thus we neglect this difference and will be working with 
correlation function given above. 



1. NUMERICAL SIMULATION. 

We simulated primordial dynamical equations ([1]) 
in a periodic spatial domain 2n x 2-7T. Main part of 
the simulations was performed on a grid consisting of 
1024 x 1024 knots. Also we performed long time simula- 
tion on the grid 256 x 256. The used numerical code [2T] 

was verified in 0[MinimHH[H[23] ' Gravit y accel " 
eration was g = 1. Pseudo- viscous damping coefficient 
had the following form 
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where kd = 256 and 70,1024 = 2.7 x 10 4 for the grid 
1024 x 1024 and k d = 64 and 70,250 = 2.4 x 10 2 for 
the smaller grid 256 x 256. Pumping was an isotropic 
driving force narrow in wavenumbers space with random 
phase: 



P% = fk^ t \h 
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here k p \ = 28, k P 2 



32 and F = 1.5 x 10~ 5 ; R % {t) 
was uniformly distributed random number in the in- 
terval (0, 2tt] for each k and t. Initial condition was 
low amplitude noise in all harmonics. Time steps were 
A£io24 = 6.7 x 10~ 4 and At 256 = 5.0 x 10~ 3 . We used 
Fourier series in the following form: 
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here N x , N y — are number of Fourier modes in x and 
y directions. 

As a results of simulation we observed [121IT5] for- 
mation of both direct and inverse cascades (Fig. [I] solid 
line). Average steepness was equal to -J (|V?7| 2 ) = 0.14. 
What is important, development of inverse cascade 
spectrum was arrested by discreteness of wavenumbers 
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grid in agreement with 10, 18, 24, 25]. Then large scale 
condensate started to form. The mechanism of conden- 
sate formation is the following. We have a flux of wave 
action (number of waves) from the pumping to the large 
scale region. This flux is due to nonlinear resonant in- 
teraction of waves. In |10) it was shown, that on a dis- 
crete grid of wavenumbers, which is typical for both fi- 
nite experimental wave tanks and computer simulations, 
resonance conditions are never fulfilled exactly. What 
makes it possible for them to exist is the finite width 
of the resonant curve due to nonlinear frequency shift. 
As a result, this thick resonant curve covers some knots 
of the wavenumbers grid |18) . The nonlinear frequency 
shift is proportional to the matrix element (coupling co- 
efficient) of interaction of waves, which is homogeneous 
function of the order 3 (with change of k it behaves 
as fc 3 ). Which means that we have good coverage of 
the grid knots in high wavenumbers, characteristic for 
the direct cascade, and worse and worse situation as we 
move toward the origin of the fc-plane. So at some stage 
the resonant interactions turn off due to discreteness of 
the wavenumbers grid and flux cannot propagate fur- 
ther. At the same time, new wave are still brought to 
this threshold scale from the pumping region. As a re- 
sult we have an accumulation of waves at some large 
scale, i. e. condensation. 

As one can see, value of wave action \a,k\ 2 at the con- 
densate region is more than order of magnitude larger 
than for closest harmonic of inverse cascade. Dynamics 
of large scales became extremely slow after this point. 
We managed to achieve downshift of condensate peak 
for one step of wavenumbers grid during long time sim- 
ulation on a small grid 256 x 256 (Fig. [TJ line with long 
dashes). As one can see we observed elongation of in- 
verse cascade interval without noticeable change of the 
slope. Unfortunately, inertial interval for inverse cas- 
cade is too short to exclude possible influence of pump- 
ing and condensate, we can roughly estimate its slope 
as fc -3 5 [T2], which is slightly less than the prediction of 
the theory of weak turbulence fc~ 23 / 6 ~ fc~ 3 - 83 . Resent 
findings |14j with significantly longer inertial interval for 
the inverse cascade, support this result and propose the 
slope fc -315 . 

Let us discuss possible reason for this deviation from 
the theory of weak turbulence. The direct cascade of en- 
ergy and inverse cascade of wave action correspond two 
Kolmogorov-type solutions of the Hasselmann kinetic 
equation for waves |20j . These solutions were derived by 
Zakharov in the middle of sixties under a few reasonable 
assumptions: media is isotropic with respect to rota- 
tions, dispersion relation is a power- like function uik k a , 
matrix element of nonlinear interaction (nonlinear cou- 
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Fig. [T] Spectra (|afc| 2 ). With condensate on the 
1024 x 1024 grid (solid); on the 256 x 256 grid with 
more developed condensate (long dashes). 

pling coefficient) for waves is a homogeneous function 
T(eki,ek2,ek3,ek4) = e^T(ki, £2, ^3, kij- Under these 
assumptions Zakharov [1 , 26 28 obtained Kolmogorov- 
like solutions corresponding to fluxes of two integrals of 
motion (energy and wave action or number of waves): 

nf = C x P^k-^-\ nf = C 2 Q^k- 2J ^- d . (5) 

Here d is a spatial dimension (d = 2 in our case). In 
the case of gravity waves on a deep water uj = \fgk 
(a = 1/2) and, apparently, /3 = 3. As a result one can 
get: 

„W = C 1 P 1 /3 fc -4 ! n (2) = C 2 Ql/3 fc -23/6, (6) 

The first spectrum corresponds to the direct cas- 
cade. In this work we are more interested in the second 

(2) 

spectrum n, describes to inverse cascade, correspond- 
ing to flux of number of waves (or wave action) from 
small scales (pumping) to larger scales. As one can see, 
inverse cascade spectrum formula directly depends on 
a, the power of the dispersion relation. So it would be 
helpful for understanding of the situation to measure 
directly the dispersion relation of waves. How one can 
do this measurement? 

In the case of linear plane wave, for harmonic 
will have just rotation of phase with constant amplitude 
a%(t) = exp(— iojfci), with circular frequency corre- 
sponding the linear dispersion relation uj{k) = \fgk~. So 
if one would write down the a^(t) as a function on time 
and then calculate Fourier transform from time domain 
to the frequency domain, the result of such transforma- 
tion will be (5-function in w-space. 

In the case of weak nonlinearity this rotation is the 
fastest process, but amplitude A^(t) will be already slow 
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function on time due to weakly nonlinear interaction of 
waves. Let us investigate how the dispersion relation 
is influenced by the nonlinear interactions in the sys- 
tem. For this we record a^(t) for different values of k 
and then calculate a Fourier transform on time, which 
leads to the function ar(u}). Because we have isotropic 
with respect to angle situation, we can limit ourselves 
with harmonics on any ray starting at the origin of the 
fc-plane. In our numerical experiment we recorded ev- 
ery other harmonics on the fc^-axis, from k = (4; 0) till 
k = (28; 0) and every tenth harmonic from k = (30; 0) 
till k = (180; 0) (for further wavenumbers the plot be- 
comes extremely noisy, also it is already range of ar- 
tificial dissipation) and calculated a Fourier transform. 
The resulting surface |a^(w)| is represented in Fig. [2j 
The observed sidebands are due to nonlinear interac- 
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Fig- (Color online) Surface |ag(a;)|. Normed to 
the maximum for every given value of k. Bright ver- 
ical line - pumping region. Sidebands, corresponding 
to the interaction with condensate, are clearly visible. 
Sidebands estimations: long dashes (black) - u k lp+ = 
uj k+kc + uj kc ; solid line (red) - = uj k -k c + u ka \ 

dashed-dotted line (olive) - to l ^ w+ — uj k +k a + ^>k c \ dot- 
ted line (brown) - to l ^ nn ~ — ujk-k c — ^fc c ; line with short 
dashes (blue) - linear dispersion relation ui k 
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tion of waves with condensate. Although the nonlinear 
process supposed to be weak, the fact, that condensate 
is more than an order of magnitude larger in amplitude 
(for harmonics further from the condensate difference 
is even bigger) makes them essential. If we consider 
the simplest case of three- wave nonlinear interaction of 
plane waves a^(r,t) = A^exp[i(fcr — wt)], we shall get 
the following nonlinear terms: 



a k a k c = A k n A L CX P ~i( w ( fc o) + w ( fc c) + i(feo + k c )r) 



k k a 



wavevector of some wave. The first term corresponds 
to the upper sideband, second - to the lower sideband. 
For example, let us suppose that we measured frequency 
spectrum of some wave with wavevector k, then it has to 
be k — fco + fcc for the upper sideband and k = ko~k c for 
the lower sideband. Let us recall that we have isotropic 
situation, which means that our condensate in fc-space 
is situated on a circle with radius k c (really it is an 
annulus of finite width). Which means that for every 
case we have projection on the k x axis changing from 
k c to — k c . As a result our sidebands has to be between 



up- _ 



= Uk-k c + ^fc c and ui k p = uj k+kc + uj kc for the 
upper sideband and between oj l k ow ~ = uJk-k c — UJ k c and 
UJ k° W+ = ^k+k,. + ^k c for the lower sideband. Of course 
for the complete consideration we need to take into ac- 
count the four wave interaction terms. Most probably 
this is the reason for not so perfect agreement of our 
primitive estimate and the measurements (Fig. [2] On 
this stage of the investigation we are satisfied with this 
simple qualitative explanation. 

If we consider spectrum in the vicinity of conden- 
sate peak, one can see that sidebands are of close am- 
plitudes with respect to central peak, corresponding to 
the oj(k) = \fgk dispersion relation (Fig. [3]). Which 
means that at least for inverse cascade region distortion 
of the dispersion relation have to be taken into account. 
For the direct cascade region, which is far from the con- 
densate scale, influence on the dispersion is almost neg- 
ligible, although the central line is slightly shifted by 
the nonlinear frequency shift. 
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Fig- El (Color online) Sections of the surface |ag(o;)|. 
Solid line (red): spectral line for k — (10; 0); dashed 
line (blue): spectral line for k = (20; 0). Dotted verti- 
cal lines of both colors correspond to linear dispersion 
relation ujk = \fgk- 



Here k c - is the wavevector of condensate and ko 



To summarize, we performed a direct numerical sim- 
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ulation of the isotropic turbulence of the surface gravity 
waves. Currently, the inverse cascade slope, observed 
in numerical simulations, is different from theoretical 
predictions of the wave turbulence theory. In order to 
investigate possible reasons we measured the dispersion 
relation for waves in the presence of condensate. It was 
shown that in the region of inverse cascade sidebands 
are of the same order of magnitude as the central line 
corresponding to the linear dispersion relation. It means 
that in the vicinity of condensate we have to take the 
influence of the condensate into account. One of the 
possible ways is to use the Bogolyubov transformation 
in order to calculate the augmented dispersion relation, 
in the same style as it was done for phonons in liquid 
Helium. This is problem for future work, because in our 
case situation is much more difficult, since condensate is 
located on the finite k c and coupling coefficient for grav- 
ity waves is immensely more complex. We hope that our 
current result can be one of the building blocks for the 
theory of inverse cascade in the presence of condensate. 
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V. V. Lebedev and I. V. Kolokolov during presentation 
of this work on the Scientific Council of the Landau In- 
stitute. The author was partially supported by the NSF 
grant 1131791 and grant NSh-6885. 2010.2. 

Also the author would like to thank developers of 
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